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Effect  of  Pressure  Gradient  on  First  Mode  of 
Instability  in  High-Speed  Boundary  Layers 
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and 

Yousef  H.  Zurigat1 

Mechanical  Engineering  Department 
University  of  Jordan 
Amman,  Jordan 

ABSTRACT 

The  effect  of  a  pressure  gradient  on  the  first  mode  of  instability  of  compressible  subsonic  and 
supersonic  boundary  layers  is  investigated  using  linear  stability  theory.  Formulations  are  presented 
for  nonsimilar  boundary-layer  mean  flow  and  linear  quasi-parallel  stability  problems  that  account 
for  variable  fluid  properties.  A  pressure  gradient  is  studied  that  generates  potential-flow  Mach 
number  distributions  at  the  edge  of  the  boundary  layer  of  the  form  Afc  =  ex* ,  where  c  is  a  constant 
and  x  is  the  dimensionless  stream  wise  distance.  Variations  are  calculated  for  the  maximum  growth 
rates  of  three-dimensional  first-mode  waves  with  different  edge  Mach  numbers  and  different  levels  of 
both  adverse  and  favorable  pressure  gradients.  A  favorable  pressure  gradient  is  shown  to  have  a 
stabilizing  effect  on  first-mode  waves.  However,  at  high  edge  Mach  numbers,  a  favorable  pressure 
gradient  becomes  less  effective  in  stabilizing  first-mode  waves.  The  frequencies  and  stream  wise  and 
span  wise  wave  numbers  that  correspond  to  the  maximum  growth  rates  of  first-mode  waves  decrease 
as  the  pressure  gradient  becomes  more  favorable  at  all  Mach  numbers  when  the  Reynolds  number 
R  =  1500  and  at  Mt  t  2  when  R  =  600.  Setting  the  Prandtl  number  to  unity  significantly  increases 
the  maximum  growth  rates  of  first-  and  second-mode  waves  at  high  Mach  numbers  compared  with 
setting  it  to  the  realistic  value  of  0.72. 
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Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  23681.  The  first  author’s  research  was 
supported  by  NASA  Langley  Research  Center,  Hampton,  VA,  under  contract  NASI- 19299. 
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I.  INTRODUCTION 


Previous  experimental  and  theoretical  studies  have  shown  that  an  adverse  pressure  gradient 
destabilizes  a  laminar  boundary  layer,  whereas  a  favorable  pressure  gradient  stabilizes  it.  These 
effects  should  be  understood  and  their  variations  quantified  with  the  large  parameter  space. 
Furthermore,  the  relations  between  a  pressure  gradient,  stability,  and  transition  make  quantifying 
the  effect  of  a  pressure  gradient  on  the  stability  of  boundary  layers  significant  for  predicting  and 
controlling  transition.  In  this  work,  we  study  the  effects  of  favorable  and  adverse  pressure  gradients 
on  the  stability  of  compressible  subsonic  and  supersonic  two-dimensional  (2-D)  boundary  layers.  The 
emphasis  is  on  evaluating  these  effects  on  first-mode  instability. 

The  effect  of  a  pressure  gradient  on  the  stability  of  three-dimensional  (3-D)  flows  is  different 
from  its  effect  on  2-D  flows.  The  presence  of  a  pressure  gradient  with  a  body  sweep  is  a  necessary 
condition  to  induce  crossflow  and,  therefore,  crossflow  instability  that  contributes  to  the  transition  to 
turbulence.  Although  modern  airplane  designs  use  swept  wings  that  induce  3-D  flows,  the  effect  of 
the  pressure  gradient  on  the  stability  of  2-D  flows  remains  important.  For  example,  near  the  wing 
box,  little  crossflow  is  present  and  the  flow  is  mostly  2-D.  In  addition,  the  flow  at  the  fuselages  and 
nacelles  at  low  angles  of  attack  is  axisymmetric  with  little  crossflow. 

Theoretical  studies  of  the  effect  of  a  pressure  gradient  on  the  stability  of  2-D  incompressible 
flows  have  considered  both  self-similar  (Falkner-Skan)  and  nonsimilar  boundary  layers.  These 
studies  showed  that  favorable  pressure  gradients  are  stabilizing  and  adverse  pressure  gradients  are 
destabilizing,  which  has  been  confirmed  by  numerous  experimental  investigations.  The  higher 
modes  of  instability1  in  high  Mach  number,  compressible  flows  makes  the  stability  of  these  flows 
more  complicated  than  that  of  low-speed  flows. 

A  favorable  pressure  gradient  in  compressible  boundary  layers  is  associated  with  compression 
and,  therefore,  with  an  increase  in  the  Mach  number  at  the  edge  Me  of  the  boundary  layer. 
Furthermore,  a  favorable  pressure  gradient  is  also  associated  with  a  convex  body.  However  an 
adverse  pressure  gradient  is  associated  with  expansion,  a  decrease  in  Me  and  a  concave  body.  The 
concavity  might  enhance  the  growth  of  Gbrtler  vortices  and,  consequently,  their  possible  interaction 
with  the  T-S  waves. 

For  compressible  flow,  only  a  few  published  studies  examine  the  effect  of  a  pressure  gradient  on 
laminar  boundary-layer  stability.  The  most  recent  studies  have  concentrated  on  the  second-mode 
instability  of  boundary  layers  with  a  pressure  gradient.  Malik2  showed  that  a  favorable  pressure 
gradient  stabilizes  the  second-mode  waves  by  reducing  the  peak  amplification  rate  and  decreasing 
the  band  of  frequencies  that  are  amplified  while  the  band  shifts  to  higher  values.  Some  of  Malik’s 
results  were  derived  from  viscous  stability  analysis  of  compressible  self-similar  boundary  layers  with 
favorable  pressure  gradients.  As  Malik  indicated,2  self-similar  solutions  in  compressible  boundary 
layers  with  a  pressure  gradient  are  not  known  to  exist  unless  the  Prandtl  number  Pr  is  set  equal  to 
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unity.  If  Pi*  is  set  equal  to  unity  in  air  boundary  layers,  the  stability  characteristics  of  the  flow  might 
change  significantly  in  comparison  with  setting  Pr  equal  to  the  realistic  value  of  0.72.  This  issue  is 
discussed  further  in  section  IV.  Malik's  stability  results  with  self-similar  boundary  layers  were 
computed  at  Mt  =  4.5  and  a  Reynolds  number  R  =  1500.  Malik  also  performed  calculations  for  a 
nonsimilar  boundary  layer  that  developed  on  the  wall  of  a  hypersonic  nozzle  with  a  suction  slot 
ahead  of  the  throat  For  that  case,  the  flow  accelerates  from  Mt  =  0.34  at  the  lip  of  the  suction  slot  to 
Mt  =  6  toward  the  exit  of  the  nozzle.  Malik  performed  a  stability  analysis  at  two  nozzle  stations 
corresponding  to  AT*  =  4.57  and  R  =  3165  and  to  Mt  =  5.45  and  R  =  3648.  In  both  cases,  a  favorable 
pressure  gradient  was  found  to  reduce  the  maximum  growth  rate  of  second-mode  waves  and  shift  the 
band  of  unstable  frequencies  toward  higher  values. 

The  effect  of  an  adverse  pressure  gradient  on  second-mode  instability  in  flow  past  a  flare  cone 
configuration  at  a  free-stream  Mach  number  M.  =  8  was  studied  by  Malik  et  al.s  The  flare  cone 
configurations  were  designed  to  produce  almost  constant  stream  wise  pressure  gradients.  The  mean 
flow  was  computed  with  both  parabolized  Navier-Stokes  and  boundary-layer  methods  and  the 
resulting  profiles  were  analyzed  with  linear  stability  theory.  The  amplification  rate  of  the  second¬ 
mode  disturbances  was  found  to  increase  with  the  adverse  pressure  gradient. 

Zurigat  et  al.4  showed  that  second-mode  waves  are  stabilized  by  favorable  pressure  gradients 
and  become  destabilized  by  adverse  pressure  gradients.  These  results  agree  with  those  of  Malik.2 
The  peak  amplification  rate  for  favorable  pressure  gradients  was  found  to  shift  toward  a  higher 
frequency,  which  also  agrees  with  the  conclusions  of  Malik2  at  Me  =  4.5, 4.57,  and  5.45.  However,  for 
flow  at  Me  =  6.5,  the  study  of  Zurigat  et  al.4  showed  that  the  peak  amplification  rate  shifts  toward  a 
lower  frequency  when  a  favorable  pressure  gradient  is  applied. 

Vignau®  (as  reported  in  Arnal  et  al.8)  analyzed  the  effect  of  adverse  pressure  gradients  for  an 
Me  that  varies  linearly  with  x.  He  considered  flows  with  M,  =  5,  5.8,  and  7.  Vignau  found  that 
adverse  pressure  gradients  destabilize  second-mode  waves  and  that  highly  adverse  pressure  gradi¬ 
ents  destabilize  the  usually  stable  third-mode  waves  (in  zero  pressure-gradient  flow  at  finite 
Reynolds  numbers).  The  experimental  study  of  Kimmel7  for  flow  past  a  fore  cone  over  a  flared  or  an 
ogive  aft  body  showed  that  at  Me  =  8,  favorable  pressure  gradients  delayed  transition  and  adverse 
pressure  gradients  promoted  transition.  At  Me  -  8  the  instability  was  dominated  by  second-mode 
waves. 

In  a  zero  pressure-gradient,  adiabatic  compressible  boundary  layer,  when  Af„  exceeds  a  certain 
value,  the  2-D  first-  and  second-mode  waves  merge  while  they  are  amplified.  Zurigat  et  al.4  found 
that  both  favorable  and  adverse  pressure  gradients  have  very  little  effect  on  the  growth  rates  in  the 
merging  region.  This  result  indicates  that  at  high  Mach  numbers,  a  very  strong  favorable  pressure 
gradient  is  necessary  to  divide  the  single  unstable  region  in  the  growth-rate  frequency  domain  into 
two  separate  regions. 
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The  study  of  Vignau®  reported  that  first-mode  waves  are  stabilized  by  adverse  pressure  gradi¬ 
ents.  This  result  is  discussed  in  section  IV.  Zurigat  et  al.4  considered  the  effect  of  pressure  gradient 
on  2-D  and  3-D  first-mode  waves  at  a  single  span  wise  wave  number  and  at  the  single  value  of  Mt  = 
2.  The  flow  and  disturbance  conditions  considered  by  Zurigat  et  al.  do  not  correspond  to  the  most 
amplified  waves.  The  present  study  considers  the  effect  of  pressure  gradient  on  first-mode  waves 
under  the  most  amplified  conditions.  We  point  out  here  that,  although  the  most  unstable  second- 
mode  waves  are  2-D,  the  most  unstable  first-mode  waves  in  supersonic  boundary  layers  are  3-D.1 


H.  FORMULATION 

We  consider  the  Navier-Stokes  equations  describing  a  3-D,  unsteady,  compressible  perfect  gas 
flowing  over  a  surface.  We  choose  a  body-oriented  orthogonal  coordinate  system  with  the  axes  x*,  y* , 
and  z*  such  that  x  is  measured  from  the  attachment  line  along  the  surface,  y*  is  measured  from  the 
surface  in  the  normal  direction,  and  z*  extends  in  the  spanwise  direction  such  that  the  coordinate 
system  is  right  handed.  We  assume  that  no  body  forces  exist  and  that  no  heat  is  aaded  (from  within 
the  flow  or  from  radiation).  The  dimensional  Navier-Stokes  equations  describing  this  flow  are  the 
continuity,  momentum,  energy,  and  state  equations 

|£-  +  V*(pV)  =  0  (1) 


*( 3V  _-?•  •  -•  z* 

p  — r-+V  -V  V  =  -V  p  +V  t 

V*  ) 

•  Dh  Dp  (  _• 


Dt  Dt 

P=p'R'T' 

where 

4=^e**^+2p*e* 
<P*  =  A*(e*  j2  +2fi*e*je*j 


.  1 

^'2 


'a* 

dt)  ax’  J 
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V*  -iu*  +jv*  +kw 
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In  the  above  equations,  r*  is  the  stress  tensor,  e*  the  rate  of  strain  tensor,  <Z>*  is  the  dissipation 
function,  and  8^  is  the  Kronecker  delta.  Furthermore,  the  indices  in  the  index  notation  have  the 
values  1,  2,  and  3  and  v\=u  ,  =  v*>  v\=w*t  *i=x*,  a£=y\  and  x£=z*.  In  equation  (3) 

D/Dt,  is  the  material  derivation  given  (for  any  quantity  s*)  by 


•V  8 


(6) 


and  in  the  index  notation 


*1 

d*i 


(7) 


for  any  quantity  s’.  In  the  above  equations,  p*  is  the  density,  t  is  time,  V*  is  the  velocity  vector 
with  its  components  u  ,  v* ,  and  w*  in  the  x* ,  y*,  and  z  directions,  respectively.  The  pressure  is 
p*,  the  enthalpy  is  h* ,  the  thermal  conductivity  is  k*  ,  the  absolute  temperature  is  T* ,  the  dynamic 
viscosity  is  p*,  the  second  coefficient  of  viscosity  is  A* ,  and  the  universal  gas  constant  is  R* .  The 
unit  vectors  i,  j,  and  k  are  in  the  x* ,  y* ,  and  z  directions,  respectively. 

For  a  perfect  gas. 


f=Jc;dr- 


(8) 


We  render  the  Navier-Stokes  equations  nondimensional  by  defining 
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II 
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where  L*  is  a  reference  length  and  the  subscript  r  denotes  a  reference  quantity.  Using  and 
expanding  equations  (8M9i)  in  equations  (1H4),  we  get 
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We  can  identify  iler,  Jfr,  and  PrT  based  on  the  reference  quantities.  In  equation  (17b),  y  =  Cp/C„  is 
the  ratio  of  specific  heats,  which  is  assumed  to  be  constant. 

A.  Mean  Flow 

We  consider  2-D  steady  flow  over  a  body  extending  infinitely  in  the  spanwise  direction.  The 
body  can  be  curved;  therefore,  it  could  exhibit  a  pressure  gradient.  However,  the  flow  is 
perpendicular  to  the  body  at  the  attachment  line  (no  sweep,  spanwise  velocity  w  =  0);  therefore,  the 
flow  is  2-D.  It  follows  that 

4p0,  ^  =  0>  “’  =  0  (18) 

Our  interest  here  is  in  boundary-layer  flows.  To  derive  the  boundary-layer  equations  for  the  flow  we 
are  considering,  we  let 


(19b) 


Furthermore,  we  take  the  reference  quantities  as  those  at  the  free  stream  °°  and  replace  the 
subscript  r  with  <*>.  Substituting  equations  (18)-(19b)  into  equations  (10)-{15)  and  assuming 
Re,  -  Rem  is  large,  we  obtain,  to  the  leading  order, 


!<*•)♦!(*)- o 

(20) 

J  du  dp  d  f 

(21) 

#=0 

(22) 

(23) 

and 

P  =  fMjT  (24) 

The  boundaiy  conditions  at  the  wall  are 

u  =  0  at  y-0  (25) 

v  =  ^Re^vu  =  jfr  at  y  =  0  (26) 

Up* 

and 


or 


T* 

T  =  TW=  — t- 


at  y  =  0 


(27) 


Equation  (25)  is  the  no-slip  condition.  In  equation  (26),  v*w  is  the  suction  (vw  <  0)  or  blowing  ( vw  >  0) 
velocity  at  the  wall  and  for  an  impermeable  wall  v  =  u*  =  0.  The  first  condition  in  equation  (27)  is  for 
an  adiabatic  surface,  whereas  the  second  is  for  some  specified  distribution  of  the  wall  absolute 
temperature  T*.  At  a  large  y  we  have 

u-*ue  as  y  (28) 

T  ~*Tt  as  y  -*<*>  (29) 
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where  u,  and  Tt  are  the  nondimensional  stream  wise  velocity  and  temperature  at  the  edge  of  the 
boundary  layer,  respectively.  They  are  made  nondimensional  with  respect  to  f/1  and  T*  and,  thus, 


It  follows  from  expression  (22)  that 
P  =  P(*) 

If  equation  (24)  is  applied  the  edge  of  the  boundary  layer,  then  we  get 

Pe  =  PJLTt 


(30a) 

(30b) 


(31) 


(32) 


At  the  same  x,p  =  pe;  therefore. 


pT  =  peTe 


(33) 


To  transform  the  boundary-layer  equations  using  the  Levy-Lees  variables,  we  let 
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With  equations  (33)  and  (34),  we  can  show  that  equations  (20H24)  transform  into 

24F{+V,  +  F  =  0 

(35) 

2 &Ft  +  VFn  -(&„)"  +  Po(F2  -Q)= 0 

(36) 

(37) 

To  reach  equation  (36),  we  used  the  relation 

dp  du. 

it=-p'u'-£ 

(38) 

which  follows  from  applying  equation  (21)  at  the  edge  of  the  boundary  layer. 

The  boundary  conditions  in  expressions  (25H29)  transform  into 

F  =  0  at  tj  — 

0 

(39) 

at  r,  =  0 

ptptU'  v  U„ 

(40) 

Q„  =  0 

or 

Q=<L  =  %- 

at  r/  =  0 

(41) 

F->  X  Q-*  1 

as  tj-^oo 

(42) 

The  variations  of  the  fluid  properties  with  temperature  are  given  in  appendix  A.  The  potential  flow 
quantities  and  their  relation  to  the  free-stream  quantities  are  given  in  appendix  B. 

For  a  constant  Cp  and  a  constant  Pr,  Cp  e  =  Cp-  1  and,  therefore,  equations  (35M37)  reduce  to 

the  familiar  form 

2^+V,  +  F  =  0 

(43) 

2 &Ft  +VFr)-(6F11)rj  +  p0(F2 -Q)  =  0 

(44) 

Z&Qs+VQ,  ~(6Q,)n  ~(7~  =0 

(45) 

The  boundary  conditions  corresponding  to  those  given  by  equations  (39H42)  do  not  change  with  a 
constant  Cp  and  Pr. 
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To  study  the  effect  of  pressure  gradients  on  the  stability  of  compressible  boundary  layers,  we 
could  select  a  certain  body  with  a  surface  that  exhibits  a  pressure  gradient,  solve  the  mean-flow 
problem  over  the  surface,  and  study  the  stability  of  the  resulting  profiles.  The  second  option  is  to 
specify  some  streamwise  variation  of  the  potential  flow  velocity  ut  or  the  potential  flow  Me.  In  the 
study  of  Zurigat  et  al.,4  as  well  as  in  this  work,  the  pressure  gradient  is  assumed  to  correspond  to  a 
boundary-layer  edge  Mach  number  distribution  of  the  form 

Me=cxn 

where  c  is  a  constant  and  x  is  the  dimensionless  streamwise  distance. 


B.  Stability 

We  consider  the  linear  quasi-parallel  stability  of  the  mean  flow  described  in  the  previous  section 
III.A  in  relation  to  3-D  unsteady  disturbances.  Stability  studies  that  account  for  the  nonparallel 
effects8-14  in  compressible  flow  over  a  flat  plate  showed  that  these  effects  can  be  significant  within 
certain  ranges  of  the  parameter  space.  The  effect  of  nonparallelism  on  the  stability  of  2-D 
compressible  flow  with  a  pressure  gradient  is  unknown;  therefore,  the  results  presented  in  this  work 
using  the  quasi-parallel  assumption  should  be  interpreted  accordingly.  Furthermore,  the  effect  of 
the  body  curvature  on  the  stability  of  2-D  high-speed  flows  is  unknown.  Preliminary  calculations 
performed  by  Malik  (private  communications,  1993)  indicate  that  this  effect  can  be  significant  for 
certain  parameters.  However,  this  effect  is  not  accounted  for  in  the  present  study.  The  stability 
analysis  accounts  for  variable  fluid  properties;  such  a  formulation  with  these  properties  was 
introduced  earlier, 1S-9'14  and  we  present  it  here  for  completeness.  To  perform  the  analysis,  we  scale 
the  Navier-Stokes  equations  and  make  them  nondimensional  by  defining 
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°r 


(46a) 


(46b) 


(46c) 

(46d) 


(46e) 
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(M= 


(iV) 


K 

K  =  ~ 

*e 


r* 

T  =  - t 

r.* 


where 


5r*=te-  (47) 

V  “« 

We  use  equations  (46a)-(46i)  in  equations  (l)-(4)  and  assume  each  dependent  quantity  as  a 
superposition  of  a  mean  quantity  and  a  disturbance  quantity;  thus,  we  let 


{u,v,w)  =  (um,um,0)  +  {u,v,w) 


P  =  Pm+P 


P  =  Pm+P 


T  =  Tm+T 


K  =  Km+K 


Cp  - Cp>m  +Cp 


If  we  use  equations  (48aH48g)  in  the  scaled  nondimensional  Navier-Stokes  equations,  subtract  the 
basic  state,  linearize  and  impose  the  quasi-parallel  assumption,  then  we  obtain  the  following 
disturbance  equations  (after  suppressing  the  tilde  for  convenience): 

dp  du  dp  dv  dpm  dw  n 
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and 


Because  ju  and  k  are  functions  of  temperature  only,  we  have 


I u  = 


_dpm 
dTm 


dK„ 


dTm 

We  seek  solutions  for  equations  (49H55)  in  the  form 

[u,v,p,T,w\  =  [Ci(y).  Ca(y).  ^(y).  Cs  {y)>  f7(:y)]exp|^  j  a  dx  +  fe-  at  Jj 
If  we  substitute  equation  (58)  into  equations  (49H55),  we  obtain 

D£3  +  -^rXf4 =0  , 


i(aun  - ftjJCx  +  Z3Dum  +  iaTm^  - ^{-pm(ra2  +  £2)£1-  aPnm(m  +  1)£7  +  i(m  +  l)apmDC3 

+DTmp'mD(l  +  ioiDTM s  +  &  +  H'mDum;6  ]  =  0  , 

»K  -  ^[*{m  +  +  imaDTM,  -  (a2  +  /J2)^3 

+rDTmp'mD£3  +  irnpDTmp'mCi  +  iap'mDum< £6  +  rpmD%  +  i(m  +  l)^»Z)f7]  =  0  , 

t(aam  -®)f7  +iPDTmp'nt;3  +  ipm{m  +  l)PD£3 

-fim(a2  +  rp2)C^DTmp'mD^+pmD%}  =  0 
and 

*>(«*  -  m)?6  +  DTmCpC3  -  i(r«  -  l)Me2rm(aum  -  co^ 

_  (Zi  [2Duniim{DS \  +  iai3 )  +  p’n{Dum  fc] 

-^rK(a2 +P%+KmD2Z*+2 K’mD\- +(DTm)2K"m;5}=0  , 
where  the  prime  indicates  differentiation  with  respect  to  Tm  and  D  =  d/dy. 
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(56d) 

(57a) 

(57b) 

(58) 

(59) 

(60) 

(61) 

(62) 

(63) 


Equations  (59H63)  can  be  cast  into  the  following  set  of  eight  coupled  first-order  ordinary  differ¬ 
ential  equations: 

8 

BC.-J/^C*^,  n  =  X2,...,8,  (64) 

m-1 

where  £2  =  i*Ci»  Ce  =  DCs>  and  C8  =D£ 7  and  the  nonzero  coefficients  are  defined  in  appendix 
C.  The  boundary  conditions  are 


Ci  -Cs  =Cs  -Cr  -0 

at  y  =  0 

(65) 

Cn  is  bounded  as 

y 

(66) 

The  corresponding  equations  for  constant  Cp  and  Pr  are  obtained  by  setting 

Cp  ~  Cp,m  -  1 

and 

yt-y-  constant 

in  the  coefficients  in  appendix  C.  The  resulting  equations  are  found  in  many  references  (e.g.,  Zurigat 
et  al.4) 

IH.  METHODS  OF  SOLUTION 

The  nonsimilar  boundary-layer  mean-flow  equations  with  pressure  gradient  are  solved  with  an 
initial  zero  pressure-gradient  region  of  length  Xq  to  avoid  the  singularity  at  x  =  0.  At  x  =  xo  (<  1),  the 
flow  starts  to  undergo  compression  or  expansion  at  a  rate  determined  by  the  choice  of  n  in  the 
equation  Me=cxn.  Positive  values  of  n  produce  favorable  pressure  gradients;  negative  values  of  n 
produce  adverse  pressure  gradients.  This  compression  or  expansion  terminates  at  x  =  I;  Me  reaches 
a  desired  value  Md  that  is  set  a  priori.  The  free-stream  Mach  number  is  chosen  such  that  Me  =  Md 
when  (x  =  1).  Stability  calculations  are  carried  out  at  x  =  1,  which  allows  a  comparison  to  be  made 
between  favorable  and  adverse  pressure  gradient  results  at  the  same  edge  Mach  number  and 
stability  Reynolds  number. 

Equations  (43M45)  and  the  boundary  conditions  in  equations  (39H42)  are  solved  with  a  finite- 
difference  method  and  a  Newton-Raphson  linearization.  The  discretization  and  solution  procedure 
follow  those  described  by  Blottner.16  The  special  case  of  self-similar  boundary  layers  is  obtained  by 
dropping  the  leading  terms  in  equations  (43H45).  This  self-similar  solution  is  then  used  to  start  its 
nonsimilar  counterpart.  The  dependence  of  the  value  of  the  pressure  gradient  parameter  on  the 
length  of  the  initial  zero  pressure-gradient  region  x0  is  very  weak.  For  example,  at  M «  =  4.5  and 
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n  =  -0.2,  a  change  of  x0  from  0.2  to  0.05  increases  the  value  of  A>  by  only  0.3  percent.  However,  a 
switch  from  a  zero  pressure  gradient  to  a  nonzero  pressure  gradient  flow  region  corresponds  to  a 
sudden  change  in  the  geometry  of  the  surface  and  could  introduce  an  upstream  effect  on  the  mean 
flow;  the  parabolic  nonsimilar  boundary-layer  theory  does  not  account  for  this  effect 

The  stability  problem  is  solved  by  a  second-order-accurate  finite-difference  scheme  with  deferred 
correction17  coupled  with  a  Newton-Raphson  iteration  on  the  eigenvalue.  More  details  of  the  method 
are  given  by  Asfar  et  al.18 

IV.  RESULTS 

A.  Mean  Flow 

Variations  of  Mt  with  the  stream  wise  distance  are  shown  in  Figure  1  for  several  values  of  n  .  As 
was  discussed  in  the  previous  section,  the  length  of  the  initial  region  in  which  the  pressure  gradient 
is  assumed  to  be  zero  has  a  very  small  effect  on  the  mean  flow  calculations  at  x  =  1.  The  stream  wise 
velocity  and  the  temperature  profiles  at  x  =  1  with  the  same  parameters  shown  in  Figure  1  are 
shown  in  Figures  2(a)  and  2(b),  respectively.  Figure  2(a)  clearly  shows  that  a  favorable  pressure 
gradient  changes  the  streamwise  velocity  profile  by  making  it  fuller,  which  stabilizes  the  flow. 

The  zero  pressure-gradient  compressible  flow  has  a  single  generalized  inflection  point  inside  the 
boundary  layer;  a  favorable  pressure  gradient  creates  another  generalized  inflection  point  near  the 
wall  (Fig.  3).  The  generalized  inflection  point  is  defined  as  the  point  where  D(pmDum )  =  0  with 
D  =  d/dy,  um  is  the  streamwise  mean-flow  velocity  component,  and  pm  is  the  mean-flow  density. 
The  significance  of  a  generalized  inflection  point  at  y  =  y,  inside  the  boundary  layer,  where  yt  >  ya 
at  which  um(y0)  =  l-M~1,  is  that  it  makes  the  boundary  layer  unstable  to  inviscid  disturbances. 
The  generalized  inflection  point  can  be  thought  of  as  the  compressible  analog  of  the  inflection  point 
in  incompressible  flows,  which  is  responsible  for  making  the  flow  inviscidly  unstable.  As  the 
favorable  pressure  gradient  level  increases,  the  two  generalized  inflection  points  move  closer  until 
they  meet  and  disappear  (Fig.  3).  Figure  3  clearly  demonstrates  that  the  minimum  favorable 
pressure  gradient  level  needed  to  eliminate  the  generalized  inflection  points  increases  as  Me 
increases. 

B.  Stability 

The  quantities  R,  F,  and  B  will  be  used  to  present  the  stability  results.  The  Reynolds  number  R 
is  given  by  equation  (56c).  The  frequency  f*  is  the  dimensional  frequency  of  the  disturbance  in 
cycles  per  second  (hertz),  and  it  remains  constant  for  the  same  physical  wave.  The  dimensional 
circular  frequency  cd *  is  given  by 

(0*  =  2ttf*  (67) 
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and  the  nondimensional  circular  frequency  a>  is  defined  as 
a>  =  a>'8'r/U' 

where  8r  is  defined  in  equation  (47).  The  frequency  parameter  F  is  defined  as 
F  =  oj/R 

and  equations  (67H69)  and  equation  (56c)  show  that 

r  *2 
“e 


(68) 


(69) 


(70) 


The  dimensional  spanwise  wave  number  is  p* ,  which  remains  constant  for  the  same  physical  wave. 
The  nondimensional  spanwise  wave  number  /)  is  defined  as 

P  =  (71) 

We  also  define  a  nondimensional  spanwise  wave  number  parameter  B  as 

P  =  1000 P  /  R  (72) 

Equations  (71),  (72),  and  (56c)  lead  to 


B  1000 P'vl 

Furthermore,  Rex  is  defined  as 


(73) 


or 


Rex  =-£t- 


or 


Rez  = 


/  •  \ 
ifs. 


Re=R2 


(74) 


We  start  by  presenting  stability  results  that  demonstrate  the  importance  of  the  first  mode  in  the 
transition  process.  To  this  end,  the  predicted  transition  location  in  compressible  flow  over  a  flat 
plate  was  calculated  from  the  empirical  e9  method.  A  wave  with  a  fixed  dimensional  frequency  and  a 
fixed  dimensional  spanwise  wave  number  was  followed  as  it  propagated  downstream.  The  Reynolds 
number  where  the  N  factor  reached  9  ((Rex)jy=9)  was  determined.  The  dimensional  frequency  and 
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spanwise  wave  number  were  varied  to  compute  the  lowest  (ReI)w„9.  Variation  of  the  predicted 
transition  Reynolds  number  with  the  free-stream  Mach  number  when  transition  is  due  to  first-  or 
second-mode  waves  is  shown  in  Figure  4.  The  figure  clearly  shows  that  the  e9  method  predicts  that 
transition  is  due  to  first-mode  waves  up  to  M_  =  6  to  6.5.  Results  similar  to  those  shown  in  Figure  4 
for  flow  past  a  cone  with  zero  and  small  pressure  gradients  were  computed  and  presented  by 
Lysenko.19  The  computations  were  performed  by  solution  of  the  Dunn-Lin20  stability  equations  and 
the  results  show  that  the  predicted  transition  is  due  to  first-mode  waves  up  to  M  *  7. 

To  study  the  effects  of  favorable  and  adverse  pressure  gradients  on  the  first  mode  of  instability 
in  compressible  boundary  layers,  we  considered  flows  with  a  pressure  gradient  parameter  n  =  0.1,  0, 
and  -0.1  and  a  potential  flow  Mt  that  ranges  from  0  to  the  hypersonic  value  of  7.4.  The  variation  of 
the  maximum  growth  rate  of  first-mode  waves  with  Me  was  computed  for  all  three  values  of  n.  The 
results  are  shown  in  Figure  5(a).  The  growth  rates  were  maximized  over  both  the  frequency  and  the 
spanwise  wave  number.  The  step  in  F  was  chosen  to  be  0.1  x  MT6,  and  the  step  in  the  spanwise 
wave  number  0  was  chosen  to  be  0.001.  The  stabilizing  effect  of  a  favorable  pressure  gradient  on  the 
first  mode  of  instability  and  at  all  values  Me  is  very  clear  in  Figure  5(a).  The  effectiveness  of  a 
favorable  pressure  gradient  in  reducing  the  maximum  growth  rates  of  first-mode  waves  diminishes 
at  high  supersonic  and  hypersonic  Mach  numbers,  which  is  also  evident  in  Figure  5(a).  Note  in 
Figure  5(a)  that  for  zero  and  adverse  pressure  gradients,  compressibility  reduces  the  maximum 
growth  rate;  for  the  favorable  pressure  gradient  case  of  n  =  0.1,  compressibility  increases  the 
maximum  growth  rate.  The  increase  in  the  maximum  growth  rate  when  n  =  0.1  results  from  the 
ineffective  pressure  gradient  at  high  values  of  Mt.  A  decrease  in  the  effective  pressure  gradient 
when  the  value  of  the  critical  Reynolds  number  changes  as  the  Mach  number  increases  was 
demonstrated  in  the  asymptotic  calculations  of  Shapiro21  for  2-D  disturbances  at  Me  =  1.5  to  3.  Also 
shown  in  figure  5(a)  is  the  variation  of  the  maximum  growth  rate  of  second-mode  waves  at  the  same 
R  and  for  the  same  values  of  it. 

The  streamwise  wave  numbers,  frequencies,  and  spanwise  wave  numbers  that  correspond  to  the 
maximum  growth  rates  in  Figure  5(a)  are  shown  in  Figures  5(b)  through  5(d),  respectively.  By 
increasing  Me,  the  most  amplified  spanwise  wave  number  remains  at  zero  (2-D  wave)  to  nearly  sonic 
values,  where  it  increases  sharply,  reaches  a  peak,  then  starts  to  decrease  slowly  (Fig.  5(d)).  Then  it 
encounters  another  sharp  increase,  reaches  a  second  peak,  and  is  followed  by  a  second  range  of  slow 
decrease.  As  the  pressure  gradient  becomes  more  favorable,  both  peaks  shift  toward  higher  values  of 
Me.  Figures  5(b)  and  5(c)  show  that  the  most  amplified  frequencies  and  streamwise  wave  numbers 
of  first-mode  waves  decrease  as  the  pressure  gradient  becomes  more  favorable.  The  most  amplified 
spanwise  wave  numbers  decrease  for  all  values  of  Mt  as  the  pressure  gradient  becomes  more 
favorable. 


17 


Results  similar  to  those  shown  in  Figure  5  but  at  R  -  600  are  shown  in  Figure  6.  The  values  of 
n  are  -0.12,  -0.06,  0,  and  0.06.  The  general  features  in  Figure  6  at  R  =  600  are  similar  to  those  in 
Figure  5  at  R  -  1500.  However,  the  most  amplified  frequencies  and  stream  wise  wave  numbers  in 
figure  6  decrease  as  the  pressure  gradient  becomes  more  favorable  at  Me  2  2,  whereas  in  figure  5 
they  decrease  at  all  values  of  M(,  although  the  decrease  in  Figure  5  is  small  at  low  values  of  Mt. 
This  performance  of  the  streamwise  wave  numbers  and  the  frequencies  that  correspond  to  the 
maximum  growth  rates  of  first-mode  waves  at  R  =  600  are  similar  to  those  of  second-mode  waves 
(Figs.  5(b)  and  5(c)).  This  result  indicates  that  if  we  compute  and  plot  the  variation  of  the  maximum 
growth  rate  (over  all  spanwise  wave  numbers)  versus  the  frequency  at  fixed  Me  and  at  R  *  600,  then 
at  low  values  of  Mt  the  peak  amplification  is  expected  to  decrease  and  shift  toward  higher  values  of 
F  as  the  pressure  gradient  becomes  more  favorable.  However,  at  high  values  of  the  peak 
amplification  is  expected  to  decrease  and  shift  toward  lower  values  of  F  as  the  pressure  gradient 
becomes  more  favorable. 

The  decrease  in  the  effective  favorable  pressure  gradient  on  stabilizing  the  first-  and  second¬ 
mode  disturbances  at  high  Mach  numbers  is  disturbing  in  regard  to  the  natural  laminar  flow  of  high¬ 
speed  flow.  This  ineffectiveness  suggests  that  other  methods  should  be  explored  to  delay  transition 
and  to  improve  the  efficiency  of  high-speed  aerodynamic  surfaces. 

In  the  results  presented  thus  far,  first-mode  waves  were  found  to  become  destabilized  by 
adverse  pressure  gradients.  These  results  seem  to  disagree  with  the  results  of  Vignau6  (as  reported 
in  Amal  et  al.6).  The  results  of  Vignau  show  that  the  critical  Reynolds  number  of  first-mode  waves 
increases  with  adverse  pressure  gradients.  This  finding  also  seems  to  disagree  with  the  results  of 
the  asymptotic  calculations  of  Shapiro31  for  two-dimensional  first-mode  disturbances.  The 
calculations  of  Shapiro23  showed  clearly  that  the  critical  Reynolds  number  decreases  with  an  adverse 
pressure  gradient  for  Mt  =  1.5  to  3.  However,  the  details  of  the  work  of  Vignau6  are  not  available  in 
English.  Furthermore,  the  reference  length  scale  in  the  work  of  Vignau  is  the  displacement 
thickness  and  not  5*  (given  by  equation  (47)).  We  believe  that  what  seems  to  be  a  disagreement 
among  the  results  of  Vignau,  the  present  results,  and  the  results  of  Shapiro  is  more  likely  due  to  the 
use  of  different  reference  length  scales  in  the  two  sets  of  results. 

Zurigat  et  al.4  showed  that  the  variation  of  the  maximum  growth  rate  of  2-D  second-mode  waves 
(for  the  considered  flow  with  Mt  =  ex")  with  n  is  almost  linear.  This  linear  performance  was  found  to 
break  down  close  to  separation.  To  determine  if  the  linearity  persists  in  the  case  of  first-mode  waves, 
the  maximum  growth  rate  of  first-mode  waves  (over  all  frequencies  and  spanwise  wave  numbers) 
was  computed  as  a  function  of  n  for  Me  =  0, 2, 4,  and  6  at  R  =  1500.  The  results  are  shown  in  Figure 
7,  and  the  variation  is  clearly  almost  linear.  However,  the  linearity  breaks  down  close  to  separation 
(not  shown  in  the  figure),  and  the  maximum  growth  rate  increases  significantly.  Such  linearity  can 
be  utilized  to  compute  the  maximum  growth  rate  for  any  n  based  on  the  values  at  two  different 
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values  of  n  and  can  also  be  used  in  the  context  of  an  envelope  method  to  correlate  transition.  The 
circles  in  Figure  7  indicate  results  of  calculations  in  which  the  length  of  the  initial  zero  pressure- 
gradient  region  is  x®  =  0.1,  whereas  for  the  rest  of  the  results,  Xq  =  0.2.  Thus,  a  change  in  Xq  from 
0.2  to  0.1  has  virtually  no  effect  on  the  value  of  the  maximum  growth  rate. 

Next,  we  considered  the  transverse  distribution  of  fluctuation  amplitude.  Computations  were 
made  at  M,  =  0,  2,  and  6  for  an  adverse,  zero,  and  a  favorable  pressure  gradient  The  parameters 
correspond  to  some  points  in  Figure  7.  For  Mt  =  0,  only  the  distribution  of  streamwise  velocity 
fluctuations  is  presented,  whereas  for  Mt  =  2  and  6,  the  distributions  of  streamwise  velocity  and 
temperature  fluctuations  are  presented.  The  results  for  Me  =  0,  2,  and  6  are  shown  in  Figures  8-10, 
respectively.  In  Figure  8  and  for  n  =  0,  the  distribution  is  standard  for  2-D  disturbances  in 
incompressible  flow  over  a  flat  plate,  which  is  characterized  by  two  peaks  wherein  the  peak  close  to 
the  wall  is  larger  than  the  peak  away  from  the  wall.  A  favorable  pressure  gradient  shifts  the  lower 
peak  slightly  and  the  upper  peak  significantly  toward  the  wall.  The  adverse  pressure  gradient  with 
a  magnitude  close  to  inducing  separation  creates  the  three-peak  structure  shown  in  Figure  8.  This 
structure  is  very  similar  to  that  noted  theoretically22-23  and  experimentally24  within  a  separation 
bubble  caused  by  a  roughness  element.  The  middle  peak  is  believed  to  be  due  to  shear-layer 
instability.  Note  in  Figure  8  that  as  the  pressure  gradient  becomes  less  adverse  or  more  favorable, 
the  disturbance  decays  faster  in  the  free  stream.  The  distribution  of  the  phase  at  n  -  -0.1  (not 
shown)  exhibits  two  phase  jumps  similar  to  those  found  theoretically26  and  experimentally24  in  a 
separation  bubble  caused  by  a  roughness  element.  For  Mt  =  2,  the  distributions  of  the  magnitudes 
of  streamwise  velocity  and  temperature  fluctuations  are  shown  in  Figures  9(a)  and  9(b),  respectively. 
For  both  quantities,  the  distributions  of  the  adverse  pressure  gradient  are  characterized  by  three 
peaks,  whereas  for  n  =  0  and  n  =  0.11,  only  two  peaks  are  evidenced.  The  temperature  fluctuation 
distribution  in  Figure  9(b)  at  n  =  -0.15  is  similar  to  that  found  by  Masad  and  Nayfeh23  for  separating 
flow  over  a  roughness  element  at  Af*,  =  0.8.  For  Mt  =  6,  the  distributions  of  the  magnitudes  of 
streamwise  velocity  and  temperature  fluctuations  are  shown  in  Figures  10(a)  and  10(b),  respectively. 
Note  that  the  shifts  in  the  peaks  in  Figures  8-10  do  not  necessarily  correspond  to  the  shifts  of  the 
dimensional  physical  fluctuation.  This  situation  occurs  because  y  in  Figures  8-10  is  y*  18*  and  8*  is 
given  by  equation  (47),  which  shows  clearly  that  8*  depends  on  the  value  of  the  pressure  gradient 
and  Mach  number.  Therefore,  the  shifts  should  be  interpreted  accordingly. 

As  mentioned  in  the  introduction,  self-similar  solutions  of  compressible  boundary  layers  with  a 
pressure  gradient  exist  only  when  the  Prandtl  number  is  set  equal  to  unity.  Furthermore,  many 
investigators  use  a  Prandtl  number  of  unity  because  it  simplifies  the  mean-flow  problems.  When  the 
Prandtl  number  is  set  to  unity  and  a  pressure  gradient  is  present ,  an  analytic  expression  exists  for 
the  temperature  profile.26  This  same  expression  holds  when  the  Prandtl  number  is  unity,  the  wall  is 
adiabatic,  and  the  pressure  gradient  is  zero.26  Changing  the  Prandtl  number  from  0.72  to  1.0  causes 
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a  relatively  small  increase  in  the  maximum  growth  rate  at  subsonic  Mach  numbers  (Fig.  11); 
however,  it  increases  the  maximum  growth  rate  significantly  at  high  Mach  numbers.  The  significant 
destabilizing  effect  of  setting  the  Prandtl  number  to  unity  at  high  Mach  numbers  is  true  for  both 
first-  and  second-mode  waves.  When  the  Prandtl  number  is  changed  from  0.72  to  1.0,  the 
frequencies  and  the  streamwise  and  spanwise  wave  numbers  that  correspond  to  the  maximum 
growth  rates  of  the  first-mode  waves  increase,  particularly  at  high  Mach  numbers. 

V.  CONCLUSIONS 

Linear  stability  theory  is  used  to  study  the  effect  of  a  pressure  gradient  on  the  first  mode  of 
instability  in  compressible  subsonic  and  supersonic  boundary  layers.  Formulations  are  presented  for 
the  nonsimilar  boundary-layer  mean  flow  and  stability  problems  that  account  for  variable  fluid 
properties.  Consideration  is  given  to  the  effects  of  both  favorable  and  adverse  pressure  gradients  on 
the  stability  of  flows  at  Mach  numbers  of  0  to  7. 

Although  a  favorable  pressure  gradient  does  have  a  stabilizing  effect  on  first-mode  waves,  it 
becomes  less  effective  in  stabilizing  these  waves  at  high-edge  Mach  numbers.  Also,  both  the 
frequencies  and  the  streamwise  and  spanwise  wave  numbers  that  correspond  to  the  maximum 
growth  rates  of  first-mode  waves  decrease  as  the  pressure  gradient  becomes  more  favorable  at  all 
Mach  numbers  when  the  Reynolds  number  is  1500  and  at  Mach  numbers  that  are  higher  than 
approximately  2  when  the  Reynolds  number  is  600.  Setting  the  Prandtl  number  to  unity  does 
increase  the  maximum  growth  rates  significantly  at  high  Mach  numbers. 
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Appendix  A  -  Variations  of  Fluid  Properties  with  Temperature 


.  1.458T*,,i  x  10-8 

U  = - ; - 

T  + 110.4 


The  variations  of  u  ,  k*  ,  und  C*  with  temperature  are  given  by 

for  T*  >  110.4, 

H*  =  0.693873 x  lO^T*,  for  T*  <>  110.4, 

.  0.6325r*12  x  10-*  , 

k  = — - — - — -r,  for  T  >80.0, 

1+(215.4/T  )xl0_12/r 
k  =  0.222964  x  lO^T*,  for  T*  £  80, 


and 


C*  =  0.0686042(00  +  ajT*  +  o27’*J  +  a3T** ), 

where  Oq  =  3.657,  ax  =  -1.272  x  10"3,  a 2  =  2.955  x  10-6,  a3  =  -1.365  x  10“9,  all  the  quantities  are 
in  centimeter-gram-second  units,  and  T*  is  in  degrees  Kelvin. 
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Appendix  B  -  Potential  and  Free-Stream  Flow  Relationships 


For  steady,  2-D  inviscid  (potential)  flow  at  the  edge  of  the  boundary  layer,  equations  (10H15) 
reduce  to 


£(P.“«)+^(p«y.)  =  ° 


(  du.  du.  \  dp. 

'  du.  du.  \  dp. 

P'[cp,c“'  +  Cp,'ve  =  ^M~  Ut^t  +  V'^!fy 


P,=PeRJT, 


where 


p  _r*t: 
”  *2 
u: 


(Bl) 

(B2) 

(B3) 

(B4) 

(B5) 

(B6) 


The  reference  quantities  are  based  on  «■  and  e  . 

By  considering  the  leading  order  and  by  ignoring  the  viscous-inviscid  interaction  problem,  we 
obtain  ve  =  =  0  even  when  suction  velocity  is  present.  Therefore,  equations  (B1MB5)  become 

CK 


and 


where 


f(pA)^  =  0 

(B7) 

(B8) 

ii 

o 

(B9) 

pfiP,tut  ^ = (r - 1 ]Mlut  ^ 

(BIO) 

P,  =  P'R-T' 

(B5) 

R  R'T- 

u: 

(B6) 
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If  equation  (B8)  in  equation  (BIO)  are  used,  then  we  obtain 

X 

fc,,.dr,  =  -(r-i)M24+Cl 

0 

at  x  =  0  and  ue  =  u* /C/*  =  1,  therefore, 


o 

for  a  constant  Cp,  Cp  t  =  1,  then, 

re  =  i+lzi3f2(i-ue2) 

Equation  (B5)  shows  us  that 

Using  equation  (B8)  in  the  above  equation  yields 

However,  equation  (B12)  shows  for  constant  Cp  that 

/  4\i/2 

The  above  two  equations  demonstrate  that 

However,  for  a  constant  Cp,  we  have 


constant  C. 


constant  C. 


constant  C_ 


*a  _«_• 


(r  -  _ — 


t:c;  uz 


Then, 


(Bll) 


(B12) 
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Appendix  C 

Coefficients  for  equation  (64) 


The  nonzero  elements  are 
Aa  =  As  =  As  ~  1  > 
fn=a2+p2-idH/Hm.Tm  , 
fn=-Dnm/nm  , 

As  =  + 1 )DTm/Tm -iaDnm /(im  +  RDun I ^Tm  , 

fa 4  =  +(m  + 1 )ytM2(m  , 

fa s  =  ~«(m  +  ty>ITn~ D{n'mDum)lnm  , 

As  =  //im  , 

fai  =  -ia  , 

As  =DTm/Tm  , 
fs4  -iyeM2a)  , 
fM  =  -id)/Tm  , 
fan  =  ~ip  » 

Ax  =  -iZ<x(rDTm  !Tm  +  2Z>/im  //^,)  , 

Aa  ■  ~iX<*  , 

As  =  Jr[-«2  -P2  +ioXlnmTm  +  rD2Tm ITm  +rDtimDTm  I nmTm]  , 
U=-ixnM2[aDum -dDTm/Tm-wDnn  f^]  , 

As  =  ix[rcd)um  / Tm  +  n^aDum  //im  - nuD/v  //z^]  , 
f*e=-iXr®ITm  , 

f„^-iXrPDTn/Tm-2ixPDnm/nm  , 

As  -  ~iXP  , 

fai  =  "2 (r«  -  W  Pr«  , 
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/ea  =-RPr«  DTnCp  !TnKm  -  2ia{y,  -  l)M2  Pre  Dumnn  /  Km  , 
fe 4  =  -»'(r. "  ^)MjRPr,(caim  -  <o)/Km  , 

f<*  =  iRPre  Cp(aun  -  a >)/TmKm  -  (y,  - 1 )M2  Pr,  ^(Dum)2 /  +  a2  +  02  -  K'mD2Tm  /  rm 

f66=-2K'mDTm/Km  , 

fa a  =  -*(w  +  l)^m  /  /  Mm  . 

fa*  =(m+  ifreM?Pa)+ipR/nm  , 

/as  =  -(m  +  tyf}<u/Tm  , 

fB7  =a2  +p2 -iaR! nmTm 

and 

faa  =  -Dnm/tim  , 

where 

Mm  =dnn/dTm,  DF  =  dF/dy,  d)  =  w-oum,  *  =  [i?/Mm -*>reM2<»]~\  r  =  2+~fL’ 

rhn 


K'M{DTnf/Kn  , 


=  -izL 

Mm 
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Figure  4.  Variation  of  predicted  transition  Reynolds  number  using  e9  method  with  free-stream 
Mach  number  for  flow  over  a  flat  plate  at  =  150  K  and  Pr  =  0.72. 


Me 

(a) 


Figure  5.  (a)  Variation  of  maximum  growth  rate  with  edge  Mach  number  for  three  levels  of  pressure 
gradient  at  R  =  1500,  Pr  =  0.72,  and  T„  =  120  K;  (b)  Corresponding  streamwise  wave  numbers; 
(c)  Corresponding  frequencies;  and  (d)  Corresponding  spanwise  wave-number  parameters. 
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(a) 


Figure  6.  (a)  Variation  of  maximum  growth  rate  of  3-D  first-mode  waves  with  edge  Mach  number 
for  several  levels  of  pressure  gradient  at  R  =  600,  Pr  =  0.72,  and  T„  =  120  K;  (b)  Corresponding 
stream  wise  wave  numbers;  (c)  Corresponding  frequencies;  and  (d)  Corresponding  spanwise  wave- 
number  parameters. 
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(b) 

Figure  6.  Continued. 
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Figure  7.  Variation  of  maximum  growth  rate  of  3-D  first-mode  waves  with  pressure-gradient 
parameter  for  several  edge  Mach  numbers  at  R  =  1500,  Pr  =  0.72,  and  T„  =  120  K. 
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(b) 

Figure  9.  Concluded. 
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(a) 

Figure  10.  Transverse  distribution  of  (a)  stream  wise  velocity  and  (b)  temperature  fluctuations  at 
Mo.  =  6,  R  =  1500,  Pr  =  0.72,  and  Too  =  120  K. 
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(b) 

Figure  10.  Concluded. 
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Figure  11.  Variation  of  maximum  growth  rate  with  edge  Mach  number  for  two  vslues  of  Prandtl 
number  atf?  =  1500,  Too  =  120  K,  and  zero  pressure  gradient. 
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